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The magnon pumping effect in ferromagnetic insulators under an external time-dependent 
transverse magnetic field is theoretically studied. Generation of a magnon current is discussed 
by calculating the magnon source term in the spin continuity equation. This term represents the 
non-conservation of magnons arising from an applied transverse magnetic field. The magnon 
source term has a resonance structure as a function of the angular frequency of the transverse 
field, and this fact is useful to enhance the pumping effect. 
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1. Introduction 

Recently a new branch of physics and nanotechnology 
called spintronics has emerged and has been attach- 
ing special attention from viewpoints of the fundamental 
science and application. The aim of spintronics is the 
control of the spin as well as charge degrees of freedom 
of electrons, and therefore establishing methods for gen- 
eration and observation of a spin current is an urgent 
issue. 

A standard way to generate a spin current is the spin 
pumping "^^^ effect in ferromagnetic- normal metal junc- 
tions. There the precession of the magnetization caused 
by an external field induces a spin current pumped into 
a normal metal. This method was theoretically proposed 
by R.H.Silsbee et.al^^ and Y.Tserkovnyak et.al,^^ and 
was confirmed experimentally by S.Mizukami et.al.^^ In 
a spin Hall system, i.e. in a nonmagnetic semiconductor, 
Kato et.al^^^ reported an observation of a spin current 
by measuring optically the spin accumulation which ap- 
pears as a result of spin currents at the edge of samples 
(GaAs and InGaAs). A critical issue in the observation 
of a spin current, however, is that a spin current is not 
generally conserved and therefore measuring spin accu- 
mulation does not necessarily indicate the detection of a 
spin current, in sharp contrast to the case of charge. Non- 
conservation of spins is represented by a spin relaxation 
torque, 7^, which appears in the spin continuity equa- 
tion. For a clear interpretation of experimental results 
on a spin current, to understand the relaxation torque is 
essential. 

The spin current means a flow of the spin angular mo- 
mentum in general, and in metals conduction electrons 
carry a spin current. In insulators, there is no conduc- 
tion electrons, but there exists an other kind of carrier, 
namely, spin- waves, which are collective motions of mag- 
netic moments. Experimentally, a spin- wave spin current, 
a spin current carried by spin- waves has already been es- 
tablished as a physical quantity. Kajiwara^^^ et.al have 
shown that a spin- wave spin current in an insulator can 
be generated and detected using direct and inverse spin- 



Hall effects. They have revealed the conversion of an elec- 
tric signal into spin- waves, and its subsequent transmis- 
sion through an insulator over macroscopic distance. The 
spin- wave spin current has a novel feature; this current 
persists for much greater distance than the conduction 
electron spin current in metals, which disappears within 
very short distance (typically micrometers). For exam- 
ple in the magnetic insulator (Y3Fe50i2), the spin- wave 
decay length can be several centimetres. 

In contrast to the experimental development, theoreti- 
cal studies so far of magnon transports are not enough to 
explain the experimental result on bulk systems. Meier 
and Loss^^^ have investigated the magnon transport in 
both ferromagnetic and antiferromagnetic materials and 
found that the spin conductance is quantized in the units 
of order {gji^Y /h'm the antiferromagnetic isotropic spin- 
1/2 chains {g is the gyromagnetic ratio, /iB is the Bohr 
magneton and h is Planck constant). Wang et.al^^^ have 
investigated a spin current carried by magnons and de- 
rived a Landauer-Biittiker-type formula for spin current 
transports. They have also studied the magnon trans- 
port properties of a two-level magnon quantum dot in 
the presence of the magnon-magnon scattering and ob- 
tained the nonlinear spin current as a function of the 
magnetochemical potential. These theoretical studies of 
magnon transports are limited to mesoscopic systems. 
From the viewpoint of spintronics, the magnon trans- 
port in a macroscopic scale is an urgent and important 
subject. 

In this paper, we focus on three dimensional ferromag- 
netic insulators. The magnon source term, T^n, arising 
from a time-dependent transverse magnetic field is de- 
rived microscopically through Heisenberg's equation of 
motion. We evaluate it by using Green's function with- 
out relying on the phenomenological equation, Landau- 
Lifshitz-Gilbert equation. This is the main aim of this 
paper. The emergence of this term is in sharp contrast 
to a charge current and represents the non-conservation 
of the magnon number. 

This paper is structured as follows. In §2.1, we repre- 
sent spin variables of a ferromagnetic insulator by boson 
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creation/annihilation operators via Hoist ein-Primakoff 
transformation. We then apply a time-dependent trans- 
verse magnetic field. This magnetic field generates the 
magnon source term, which breaks the magnon conser- 
vation law in the spin continuity equation. In §2.2 and 
§2.3, by evaluating the magnon source term at the low- 
temperature limit, the dependence of the magnon source 
term on the angular frequency of a transverse magnetic 
field is calculated. The magnon source term has a res- 
onating behavior when the angular frequency of an ex- 
ternal transverse magnetic field is tuned. In §2.4, the 
temperature dependence of the magnon source term is 
argued. Through the analogy with the usual conduc- 
tion electrons' spin pumping effect, the possibility for 
the magnon pumping is discussed in §3. 

2. Magnon Source Term 

2. 1 Definition 

We consider a ferromagnetic Heisenberg model in 
three dimensions. It reduces to a free boson system via 
Holstein-Primakoff transformation if we approximate the 
spin as [S : the length of a spin) ; = S — a^a^ = 
S- = (2^)i/2at[l - ata/(25)]i/2 ^ {2Sy/^al Here op- 
erators /a are magnon creation/annihilation operators 
satisfying the bosonic commutation relation. Therefore 
in the continuous limit, a three-dimensional ferromag- 
netic insulator with an external magnetic field along the 
quantization (z) axis is described at low magnon density 
limit as 



are defined as 
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Here mmag is the effective mass of a magnon and it is 
represented by a ferromagnetic exchange coupling con- 
stant in the discrete model, J, and the (square) lattice 
constant, ao, as ^^/(2mmag) = 2JSaQ^. In eq.(l), ^ is a 
constant external magnetic field along the quantization 
axis (z-axis), g is ^-factor and /iB is Bohr magneton. 
From now on including ^-factor and Bohr magneton, we 
write an external magnetic field as gfiB^ = B. We then 
apply a time-dependent transverse magnetic field with 
an angular frequency, and a constant field strength. 
To, to X-axis as T{t) = Focosl^t. 



Vr{t)=T{t) / d^x5^(x) 



a{x.t) + a^{x.t) 



(2) 



The total Hamiltonian is H = Hq -\- Vr{t). 

The magnon density, Pni(x), of the system is defined as 
the expectation value of the number operator of magnons 

Pm(x, t) = (a''"(x, t)a(x, t)). (3) 

Through Heisenberg's equation of motion, the magnon 
current density, jm, and the magnon source term, 7^, 



dpm 
dt 



(4) 



Here the magnon current density arises from the free 
part; [pm, 'Wol/l^^) = —V • jm- It reads 



J^^(x,t) = 
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-Re 
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mag 



i < {d^a\x.t))a{x.t) > 



(5) 



where /i is a direction for a magnon current to flow 
(/i = X, z). The magnon source term, which represents 
the breaking of magnon conservation, arises from a trans- 
verse magnetic field as [pm^Vr]/{ih) = 7^, i.e.. 



%n{t) = -M!^! Ini(r(t)a(xt)). 



(6) 



From now on, we treat Vr{t) as a perturbation (i.e. a 
weak transverse magnetic field) and study the effects of a 
time-dependent transverse magnetic field to the magnon 
source term. 

2.2 Evaluation 

Through the standard procedure of the Keldysh (or 
contour-ordered) Green's function, ^^^^^ the Langreth 
method, the magnon source term is evaluated (see 
also APPENDIX ) as 



(7) 



<T{t)a{t) >= J d^x' j dt'T{t)T{t'){^f/^G'(t,t') 

+0(F^). 

Here is the retarded Green's function and we have 
neglected terms which are third-order in F, which is jus- 
tified at the low magnon density regime. 

The retarded Green's function is G^iyr'^tt') = 

{^/y) Ek /(^icj/27r)e^*^-(^-^')-^^^*-*')G^k,u;' and 
G^(k, cj) = [hw — uo\^ ^ ih / {2t)]~'^ . Here V is a volume 
of the system. The lifetime r represents the damping of 
spins (r is inversely proportional to the Gilbert damping 
parameter, a). The energy cjk corresponds to the free 
part, Ho, and therefore cjk = Dk^-\-B^ D = ^^/(2mmag)- 
Then the magnon source term is calculated as 



F^ 

^ Q 



2r 



Tf cos21]t 
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2r 



{hn - B)sin21^t + ^cos2nt 



(8) 



The time average of Tm becomes 



r2 
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2r 
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2r 



4hV [(n^] + 5)2 + (A)2 5)2^(^)2 J 



(9) 



It is clear that Tm is positive (for finite temperature, see 
eq.(16) in §2.4). 
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2.3 Resonance 

We define a dimensionless quantity 7^^^ as 

1 1 



2tB\2 



4hV 



(10) 



This siiows that the magnon source term has a resonance 
structure with a time-dependent transverse magnetic 
field when the angular frequency is tuned as ^} = B/h 
(see Fig.l). This resonance is useful for the enhancement 
of the magnon pumping. 




Fig. 1. A graph of which represents the angular frequency 

dependence of the magnon source term. Parameters we have used 
are, 2t1^ = t,T = 2x 10-6[s], B = 1[G], g = 1. Therefore 2rB/h 
is 35. This T'^^ has a sharp peak around t = 35. This fact means 
that the magnon source term has a resonance structure with the 
apphed transverse magnetic field. 



2.4 Temperature dependence 

Let us look into the temperature dependence of the 
magnon source term. To do this, we include the interac- 
tion of third-order in magnon operators. Then Vr{t) is 
rewritten as 

Vr{t) = T{t) J d^xS^^) 

.r(.)/A(f)./^{„M,.,,.., 



1 

4S 



a\xt)a^xt)a{xt) + a'^ {xt)a{xt)a{xt) |, 

(11 

and the magnon source term reads 

at(xt)at(xt)a(xt) 



Tm{t) = -^^^ lm{T{t) 



a{xt) 



4S 



(12) 



Eq.(12) is calculated as 

< r{t)a{t) > = J £x' J dt'r{t)r{t')[{^y/^G'{t,t') 



(13) 



< ^aHt)aHt)a{t) > 



cfx' / dt' 



2(25)1/2 

• r{t)r{t')G^(t', t)G<{t,t) + 0(r3). 

(14) 

Here and are the advanced and lesser Green's 
functions, respectively. Because we focus on the behav- 
ior of the magnon source term at low temperature that 
we have neglected higher terms than the fourth-order in 
respect to magnon creation/annihilation operators. The 
Fourier transform of the lesser Green's function satis- 
fies, G<(k,cj) = -/B(cJk)[G'^(k,a;) - G^(k,cj)], where 
G^(k, uo) = [fvjj — cjk — ih/ (2r)]~^. Then Tm. is calculated 
as 

1 



Tm 



f2 r 
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mation 



^^{hVt^ B)sm2VLi + ^cos21]t 
^^{hQ.- B)sm2Q.i + ^cos21]t 

{n^i-BY^{l^Y . 

is k^T (/cb : Boltzmann 
Here we have approximated the 
over the Bose distribution function 



(15) 



[(2^) VV^] Ek /B('^k) = 27r dke [e^iDk^+B) - 1]" 
The time average of Tm becomes 
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(16) 



where 

r-.5-(i + J^.-»)e-»(|g)- im 

is the temperature-dependent part and 7^^^^ is defined in 
eq.(lO). When the temperature gets higher, the magnon 
source term decreases. This means that quantum and 
thermal fluctuations act in the opposite way, namely to 
increase and decrease the magnon source term, respec- 
tively. It is clear that eq.(16) reduces at the zero temper- 
ature to eq.(9). 
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3. Magnon Pumping 

The spin pumping is an effect widely used to cre- 
ate a spin current by use of the magnetization preces- 
sion. Experiments have been carried out in junctions 
of metahic ferromagnets and nonmagnetic metals. Ac- 
cording to the theory by Y.Tserkovnyak and A.Brataas,^^ 
the spin current pumped through the junction reads 
jpump _ [h/{A7T)][A,m X [dm/dt) - Ai{dm/dt)], where 
m is the magnetization direction of a localized spin and 
Ar^A\ are the interface parameters. This result is un- 
derstood by considering the spin continuity equation, 
V • js = —{dm/dt) + 7^, where js is the spin current 
density and Ts is the spin relaxation torque. In fact, the 
spin continuity equation indicates that a spin current 
is generated when the magnetization is dynamic and/or 
when Ts is finite. It has been well-known^^^ that the 
main term of spin relaxation torque, 7^, has the form 
7^ oc m X (dm/dt) in metals. As Brataas et al.^^^ have 
pointed out, this torque Ts is equivalent to the maximum 
spin current that can be drawn from the spin battery. 
The pumping formula for Jp^^^p in metals is thus under- 
stood from the spin continuity equation. The relaxation 
torque plays an essential role in spin pumping in met- 
als, and it is expected to be dominant also in the spin 
pumping in insulating ferromagnets. Our calculation of 
the relaxation torque term (i.e. the magnon source term) 
thus describes the spin pumping effect in insulators. 

We have revealed that the magnon source term has a 
sharp peak around 1] = as a result of the resonance 
with a time-dependent transverse magnetic field when 
the angular frequency is appropriately adjusted as 1^ = 
B/h. This fact is useful to enhance the magnon pumping 
effect because the external magnetic field, and the 
angular frequency of a transverse magnetic field, 1], is 
under our control. 

Experimentally, the magnon pumping effect we have 
discussed can be identified by observing the tempera- 
ture dependence of the pumped magnon current when B 
is zero; (the pumped magnon current) (x Tm oc (Ai — 
A2T^/^), where Ai and A2 are constants. 

4. Summary 

We have studied theoretically the magnon source term 
which represents the breaking of the magnon conser- 
vation law. We have revealed that the magnon source 
term has a resonance structure with an external time- 
dependent transverse magnetic field when the angular 
frequency of the applied magnetic field is tuned. This 
fact will be useful to enhance the magnon pumping ef- 
fect in insulators. This magnon pumping effect is a new 
method for a generation of a magnon current (spin- wave 
spin current) without the gradient of an external mag- 
netic field. 
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Fig. 2. (Color online) A schematic picture of the magnon pump- 
ing by a time-dependent transverse magnetic field. Thick and 
short solid arrows represent localized spins, and thin and long 
ones external magnetic fields. Circles represent magnons. By way 
of resonance with a time-dependent transverse magnetic field, the 
magnon source term is enhanced, and an enhanced magnon cur- 
rent is pumped from the ferromagnetic insulator to the adjacent 
nonmagnetic insulator. 
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Appendix: Calculation of Equation. (7) 

In this section, we briefiy show the Langreth method, 
which is useful to evaluate the perturbation expansion of 
the Keldysh ( or contour-ordered) Green's function. For 
simplicity here, we evaluate the magnon source term at 
zero temperature, eq.(7), as an example; 



a(xt) + a^xt) 



%n{t) = -^^^ lm(r{t)a{xt) ). 



(A-1) 



We have only to estimate < a(xt) >. It is evaluated as 
< a(r) > = ^Tc a(r)exp - i j^drVr{r) 

(A-2) 



j d^x'X. 



Here Tc is the path-ordering operator defined on the 
Keldysh contour, c (see Fig.A-1). We express the Keldysh 
contour as a sum of the forward path, c^, and the back- 
ward path, c^; c = c^ + c^. The integral on the Keldysh 
contour of eq.(A-2), X, is executed by taking an identity 
into account 



/ (ir^ = / dr^ + / 



(A-3) 
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as 



(A-4) 



Here is the time-ordered Green's function. By using 
the relation, G'{t,t') = G\t,t') - G<(t,t'), we obtain 
eq.(7). 



T 



T 

-X- 



(c =c^+ c^) 



Fig. A-1. Keldysh contour, c. We have taken r on forward path, 
c^. Even when r is located on backward path, c^, the result of 
this calculation is invariant because each Green's function, G^, 
G^, G^, (the greater Green's function ), is not independent; 
they obey, C^— = G^ —G^. Both the forward and backward 
paths are actually on the real axis but shifted slightly upwards 
and downwards, respectively, to distinguish them clearly. 
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